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1 Introduction 

One of the main issues of models in equilibrium statistical mechanics is to 
explain phase transitions and critical phenomena. In some sense the macro¬ 
scopic behaviour does not depend much on the microscopic one which leads 
to universal behaviour for different systems at large scales. 

In non-equilibrium systems a greater variety of phenomena occurs and 
in some sense is more sensitive to microscopic information. A unified ther¬ 
modynamical theory is much more difficult to obtain and different attempts 
were made. 

The simplest situation where one can create a so-called non-equilibrium 
steady states (NESS) is to consider a system coupled at the boundaries in 
with different reservoirs creating currents. To construct such models, one 
usually starts from a simple bulk model to which boundary terms are added. 
By simple bulk model here we mean a reversible Markovian conservative 
dynamics (such as the simple symmetric exclusion process) having simple 
stationary measures e.g. of product nature. Upon coupling such a system 
to different reservoirs, the nature of the stationary measures changes dra¬ 
matically, i.e., NESS thus created are far from product, and generically show 
long-range correlations. In some special systems, these correlations are acces¬ 
sible analytically [7j. Such models could be considered as the non-equilibrium 
analogue of exactly solvable models in equilibrium (such as the Ising model). 
On the macro scale in a large class of one-dimensional conservative Markov 
dynamics, it is believed [1] that two-point, correlation functions in NESS are 
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multilinear. On the micro scale this multi-linearity is rare and a possible 
indication of exact solvability e.g. in the sense of matrix ansatz solution [0]. 

hi 0 a class of wealth distribution models with two agents in equilibrium 
is considered. The authors showed under which conditions on the redistribu¬ 
tion measure and transition operator there exist stationary product measures 
and duality functions. At each time step two agents are redistributing a ran¬ 
dom amount e of their wealth and a constant part A G [0,1] is retained. It 
turns out that there exist only stationary product measures for the trivial 
case A = 0 and the law of e is independent of the amount of total wealth and 
in particular is Beta(k, k) distributed. 

In [I] the authors study particle systems out of equilibrium. Multilinear¬ 
ity of the two-point functions is obtained for the KMP (Kipnis-Marchioro- 
Presutti) model of heat conduction and boundary driven SEP (symmetric 
exclusion process) for trivial A, e uniform and the reservoirs follow an ex¬ 
ponential distribution. The two-point correlation functions turn out to be 
negatively correlated for the SEP resp. positively correlated for KMP while 
the off-diagonal terms are essentially the same in the large N limit. 

An example where one does not obtain multilinearity for a model of heat 
conduction can be found in [B]. The authors study the Brownian momentum 
process which is weakly coupled to heat baths and in particular the NESS 
and its proximity to the local equilibrium measure in terms of the strength 
of coupling. For three- and four-site systems, they obtain the two-point 
correlation function and show it is generically not multilinear. 

Here we are interested in the question under which conditions we have 
multilinear two-point, functions in this class of models. The conditions will 
only include a relation between the moments of the reservoir laws, moments 
of the redistribution parameter and the constant A. The result presents a 
first example where one can obtain multilinear two-point functions even in 
the absence of product stationary measures. 

The rest of our paper is organized as follows. In section 0] we will provide 
all necessary definitions. Section [3] deals with the main result and finally 
we discuss in section 4 some generalizations. Detailed computations can 
furthermore be found in the appendix. 


2 Notation and Definitions 

We consider one-dimensional interacting models on the set {1,..., N} coupled 
to some reservoirs. The sites 1,..., N can be interpreted as particles or agents. 
The boundary sites are interacting with reservoirs which are represented by 
the ghost sites 0 and N + 1. Let = S N = [0, oo)^ denote the state space 
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of the model. Each element (xi,x^) E fl can be seen as e.g. the wealth 
or energy at vertex 1, N, see 0 15] • 

For some random e E [0,1] and constant A E (0,1) we can define a map 
T^'l : Q — >■ Q for some k,l E {1,IV} with k < l by 

x’ e {Xl, -,X N ) = 

(xi ,..., x k -i, \x k + e(l - X)(x k + Xi ),..., Xxi + (1 - e)(l - X)(x k + x t ),..., x N ) 

( 1 ) 

At each time step some random amount of energy or wealth is exchanged. 
During the exchange some fixed amount A is kept while some random amount 
(modeled by e) is exchanged. 

We remark that the map conserves the total wealth (energy) s = X\ + 

... + x N . 

Let v be the law of e G [0,1] which will be called the redistribution 
parameter. We further assume (Al) that v is symmetric and its first moment 
is equal to | and that it has a second moment (A2) which will be abbreviated 
by £ e(l — e)v(de) =: a. Note that it follows trivially that by (Al), a <\. 

After some exponential waiting time with mean 1, the initial state (aq, ...,xn) £ 
Q is is redistributed and replaced by T^ ,l e (x i,..., x^). 

The generator of the model can be expressed as the sum of the bulk 
generator describing the dynamics in the bulk and boundary generators which 
describe the interaction with the reservoirs, 


s? = & L + & h + se R . 


( 2 ) 


It is a Markov process. Let / be some bounded continuous function, /:!!—>■ 
R then 


J 2 ? b (/(xi, ...,x N )) = 


f( T \!e( x i’-,x N ))v{de) - f(x i,..., x N ) 


(3) 


where p(k, l ) is the transition probability of a symmetric nearest neighbour 
random walk given by 

p{k) 0 i T 

for k, l E {2, ...,1V— 1}. The generators of the reservoirs are defined as 


22?l(/(xi, ...,x n )) = 




/(Axi + e(l - A)(a : 0 + x 1 ),x 2 , ..,x N )u(de)n L (dx 0 ) - /(xi...., x N ) 


(4) 
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resp. 


^oo /•! 


&R(f(Xi, ...,X N )) = 

L 

f(x 1, ■■■, xjv-i, Xx N + e(l - X)(x N + x N+l ))v(de)p R (dx N+l ) (5) 


1 0 JO 


- f(x 1, ...,X N ) 

/i 0 resp. [In+i denote the distributions of the reservoirs at the ghost sites 0 
resp. TV + 1. We abbreviate their first moments by 


Tl / XofiL(dxo) , T R / XN + i[i R (dx jv+i) 
Jo Jo 

resp. the second moments by 


L 2 : = 


POO 

I xlfi L (dx 0 ) , R 2 := / x^ + 1 /iR(cfeiv + i) 

Jo Jo 

We remark that at this point we only assume that the first and second mo¬ 
ments of the reservoir measures exist. 

We call a probability measure n G &(Q) on the set of probability mea¬ 
sures on stationary if and only if for / bounded and continuous on fk 


( 6 ) 


(7) 


%{f{x i, x N ))p,{dx i,..., dx N ) = 0 . 


( 8 ) 


For i < j and i,j G {0, ...,N + 1} the two-point function /i(xi,Xj) is equal 
to n(xi,Xj) := f n XiXjfj,(dxi ,..., dxjv). Later we will use for convenience the 
notation /j, tJ . Further let 

C N (i,j) ■= ti(xi,Xj) - E N (i)E N (j ) (9) 

denote the two-point correlation function of /i. We set 


resp. for 


Cjv(0,0):= / (x 0 - T L ) 2 n 0 (dx 0 ) 

Jo 


Cn(N + 1, TV + 1) / (xat + i — T R ) 2 i_i N+ i(dx n + i). 


( 10 ) 


( 11 ) 


Note that this is equivalent to saying L 2 =: /i(x 0 , x 0 ) resp. R 2 =: /z(xjv+i, Xjv+i)- 
For 0 < i < N and 1 < j < N + 1 we set Cn{ 0, j ) := CV(i, TV + 1) := 0. 

We say that the two-point function satisfies the multilinearity ansatz if 
and only if we can find some coefficients a,b,c,d,e, f, g such that for i, j G 
{0,..., IV + 1} and i < j 


Ah ,7 = a + bi + cj + dij if i < j,i = 1,..., TV - 1, j = 2,..., TV 


= e + fi + gi 2 


if i = j,i = 1, ...,7V. 


( 12 ) 


3 Result 


We present our main theorem. 

Theorem 3.1. The two-point functions are multilinear with coeffi¬ 

cients given by 

a := T 2 l 

, C Tr - T l ) [{N + 1)(A + 2a(l - X))T l + (1 - A)(l - Aa)T R ) 

(N + l)(l + N\ + 2(N-l)(l-\)a) 

(Tr - T l )T l 
c ' N + 1 

(A + 2q(1 — A)) (Tl — Tr) 2 
' (N + 1)(1 + NX + 2(N — 1)(1 — A)a) 

(1 - 2a(l - A)) (T r - Tf) [(1 + (2 N + 1)A + 4AW(1 - A ))T L + (1 - 4a)(l - X)T R ] 
1 ' (N + 1)(1 + NX + 2(N - 1)(1 - A)a)(A + 2a(l - A)) 

(l-2a(l-A ))(T l -T r ) 2 
9 ' (N + 1)(1 + NX + 2(N — 1)(1 — A)a) ' 

(13) 


if and only if L 2 = L 2 (a, A) and R 2 = R 2 (a, A) are chosen in the following 
way 

2. (1 — 2a(l — A)) 2 


L : = 


R z : = 
5 + 


_ _ r n£ 

(A + 2a(l - A)) L 

a (l — A)(l — 2a(l — \)){T l — Tr) 2 
(N + 1)[1 + NX + 2(N - 1)(1 - A)a](A + 2a(l - A)) 

rp2 

-L E 


2 (1 2a(l A)) 2 


(14) 


(A + 2a(l - X))^ R 

a( 1 - A)(l - 2a(1 - A ))(T L 


Tr) 2 


(N + 1)[1 + ATA + 2(N - 1)(1 - A)a](A + 2a(l - A)) 
In particular the two-point correlation functions are equal to 


( (l-4a)(l-A )(T L -T R f , 

\ (l+AiV+2(iV-1)(L—A)a ) j N +1 

(l-4a)(l-A)r» 


1 - 


+ 


N +1 

a(l-2a(l-A))(l-A)(r L -T fl ) 2 


0<i<j<N+l 


C N (fj) = < 


A+2a(l-A) ^ (A+2a(l-A))[l+AJV+2(JV-l)(l-A)a](iV+l) 

(l-4a)(l-A)(T fl -T L )((l+2AJV+2a(2JV-l)(l-A))r L +(l-2a(l-A))r fl ) 

(A+2<k( 1—A))(l+AJV+2(iV — I)(1—A)o ) / N+l 


( (l-4a)(l-A)(T L -r fl )^ \ i 2 . a < i = • < ^ + 1 

(l+AJV+2(JV-l)(l-A)a) J (N+l) 2 ' U — L J — JV ^ 1 


( 15 ) 
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Proof. 

We are looking for expressions for (fJ-ij)ij for i,j E {0 ,...,N+ 1}. In 
the bulk for i < j, fx^j have to satisfy the following set of equations (see 
appendix): 


{ —2(1 — A)/in + A(/Xj_iy_i + 1 + = 0 

/h+ij + /h-ij + /hj+i + Hi,j -i — 4 Hi,j = 0 

+ /h+l,j+l) + |(/h-l,i+l + Hi,i+2 ) ~~ (1 + = 0 

with A, B, C defined in the previous section as 


if z = 2,...,N- 1 

i=2,..,N—3 

j=4,...,N-l 

if i = 2,..,N- 2 
(16) 


A 

5 

(7 


l-2a(l-A) 
— + a(l - A) 


(17) 


On the boundaries we need the following conditions to be satisfied. First for 
the left reservoir (z = 1): 


f —2(1 — A) fin + A(L 2 + fi 22 ) + B(fi 1,2 + 7l1?7v(1)) = 0 

\ h'2,j + TiEn{]) + — 4/iij = 0 

[(7(Ahi + ^ 22 ) + 1(^13 + TlEn( 2)) — (1 + B)fii2 = 0 
and second on the right for j = N, 


if 3 = 1 

if j = 3,AT — 1 
if j = 2 

(18) 


f—2(1 — A)fi NN + A^jv-giv-i + -R 2 ) + B{fx N _ i jN + T r E n (N )) = 0 if z = A r 

\ /h+l,iV + /Xi- 1 ,N + T R En(i ) + Hi t N-\ — 4/ijiV = 0 if z = 2, ./V — 

[(7(/zjv-i,at-i + Maw) + \(hn- 2 ,n + T r En(N — 1)) — (1 + B)(in-i,n = 0 if z = Af — 1 

(19) 

and 

E N {i) =T L ^1 — + Tr ^ + f ) * — 0) •••>-^ + 1- (20) 

The result follows from some tedious computations. □ 


4 Conclusion and Discussion 

In this note we obtained conditions under which two-point functions in the 
non-equilibrium case will be multilinear for some class of one-dimensional 
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interacting models of wealth (heat) transport. In particular we present a 
first example of models for which the two-point function is multilinear when 
the stationary measure is not product. 

Let us make some final remarks. 

Remark I: For A = 0 and a = | we find the two-point correlation functions 
obtained for the KMP model in [T] . 

Remark II: For a non-degenerate distribution of the redistribution parame¬ 
ter e, a < Hence the occupation variables are always positively correlated. 
If v = < 5 i /2 is a degenerate measure, a = and the correlation functions are 
0 . 

Remark III: We tried to apply the multilinearity ansatz in the case that 
at the reservoirs are additionally depending on factors 7 R and 7 #, hence the 
temperature profile is not entirely linear. The boundary generators are then 
given by 


-,x N )) = 

POO pi 

7 L / /(Axi + e(l - A)(x 0 + Xi),x 2 , ..,x N )iy(de)fi L (dx 0 ) - f(x 1 ....,x N ) 

Jo Jo 

J2? R (/(X1, ...,X N )) = 


7R 



..., Zjv-i, Xx N + e(l - X)(x N + x N+1 ))u(de)n R (dx N+1 ) 


- f(x 1 , ...,x N ) 


( 21 ) 


It turns out that the two-point, functions will never be multilinear as long as 
^ 1 or T l ^ T r . 

Remark IV: Further we tried under what conditions three-point, function 
might be multilinear. Unfortunately this problem is very complex and we 
could not find a general solution to this problem even for N = 6 . 


5 Appendix 

5.1 Generator for the two-point function 

In the following we determine the generator ££ (see (J2]) ) of the two-point 
function fij(xi ,..., x^) := XiXj. It will also depend on the 1-point linear 
functions on the borders fi(x\, ■■■,Xn) := 07 . Recall that the transition op¬ 
erator T^ ,l e was defined in (JTJ) , e satisfies (Al) and (A 2 ) and the reservoirs 
have finite first and second moments.. We will distinguish 3 cases. Case I 
represents i = j, case II \i — j\ > 1 and finally case II | i — j\ = 1. 
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5.1.1 Temperature profile 

We will determine the closed form of the linear functions corresponding to 
the temperature profile of the system. Let /i be a stationary measure. We 
calculate the density profile E^{i) '■= J n Xifj,(dx i,... ,dxN ) by solving 


^f(xj)/i(dxi,..., dx jv) = 0 


( 22 ) 


for % = 1 , Let i = 2,N — 1 then we need to solve 

/ J§?(xi)/i(otei,..., dx N ) = E N (i + 1) + E N (i - 1) - 2 E N (i)) = 0 


for i = 1 and i = N we have 

[ Jf(xi)n(dxi ,..., dxAr) = 1 X (E n ( 2) + T l - 2 E n (1)) = 0 
do ^ 

/ &(x N )n(dx lt .... <fa,v) = 1—2 (T r + E n (N - 1) - 2Ejv(i)) = 0. 


(23) 


We can compute explicitely the closed form expression namely 


En{i) —T l ( 1 — 




R- 


, i — 0,..., iV -f- 1. 


(24) 


jv + iy jv + i 

Note that the linear profile is the same as in [T] for the SEP and KMP models. 


5.1.2 Case I: i—j 

We will determine the generator 

^(fii(x)) = 22? L (/ n (x)) + Sf b (fu(x)) + ^rUnnW) 

in the case that i = j. Let us first fix i = 2,...,1V — 1. First we calculate the 
generator in the bulk S^h{f-n{x))- 


-&b(fu(x)) = P( k -> 1 ) [ f( T x!e( x i^-y x N))^(de)-f(x 1 ,...,x N ) 

= T 2 (fdT” M (x)) - U(x)) + i ~ /«(*)) 

+ i(/«(r‘ +u W) - /«(*)) + i(/«(T i -M(*)) - f„(x)) 

= / ii (T i ’'+ 1 M) + MX^Hx)) - 2/ii(x) 


( 25 ) 









The boundary generators are obtained in the following way, let i— i\ 


&L(fn(x)) = fn(T 0,1 (x)) - f n (x) + fi 2 (T lt2 {x)) - fn(x) 

where 

/»oo 

fn{T 0,1 {x)) = / (Axi + e(l - A)(xi + x 0 ))v(de)n L (dx 0 ) 

Jo Jo 

= Xxi + 1(1 - A)(xi +T l ) ( 26 ) 

= + |(1 - A )T L 

and analogously for the right boundary generator S£r with Tl replaced by 
Tr and x\ by xn- Hence for i = 1 ,N the generator 2z? acting on f n (x) is 
equal to 


^(fn(x)) = J? L (fu(x)) + S?b(fii(x)) + %’n(fNN{x )) 


with 


^l(/u(*)) = (1-A) 

-W«(aO) = 


-2(1 — A)fi t i(x) + A(L 2 + f 2 , 2 (x)) + B(fi y2 (x) + T L fi(x)) 


(1-A) 


-2(1 - A)fu(x) + A(/ i _i )i _i(a;) + f i+1 , i+l (x)) + B(f i>i+1 (x) + 


%R.(fNN{x )) — 


(1-A) 


-2(1 — A)fm t ]y(x) + A(fN-i,N-i(x) + R 2 ) + -B(/jv- i,at(^) + Tr/at(x)) 


(27) 


with coefficients given by 


A = A(a, A) = - aj (1 - A) 

B = B(a, A) = A + (1 — 2a)(1 — A) 


(28) 


5.1.3 Case II: |i — j\ > 1 

We calculate 2z?(/jj(x)) for fij(x) = x^Xy Let us first assume i < j, i = 
2 ,N — 2 and j — 4,..., N — 1. We can easily see doing similar calculations 
as in the first case that 
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-%(/,iW) = / ij (r' i+1 (i))+/ iJ (r- 1 '‘(x))+/,: J .(T« +1 (x)) + /, J (T>- 1 ^w)-4/^w 

(29) 


At the left reservoir i — 1.2 < j we have 

■%(/i J (^)) = /u(r 1 ' 2 (x))+/u(T 0 ' 1 (x)) + / u (TJ9+ 1 (x)) + / l3 .(TJ- 1 9(x))-4/ u (x) 

and for the right boundary j = N, i < N — 1 


& R (fiAz)) = /*7v(r’* +1 (x))+/^(r- 1 ’ i (x))+/^(T^ +1 (x))+/ iiV (r 7V - 1 ’ Af (x))-4/ l7 v(x) 

ft follows for i < j, i = 2 ,N — 2 and j = 4,iV — 1 that the generator 
££ acting on fij(x) with \i — j\ > 1 can be written as 


•&(fij( x )) = ■Z’dfijix)) + Jf b (f i:j (x)) + Sf R {f iN {x)) 


& L {fij{x)) 

^b(Mx)) 

&R(fiN(x)) 


(1-A) 

2 

(1-A) 

2 

(1-A) 


f 2 ,j(x) + T L fj(x) + fi >j+ i(x) + fi :j -i(x) - 4 fij{x) 

fi+ij(x) + fi-i,j(x) + fi,j+i(x) + fi,j-i(x) - 4 fij(x) 

fi+l, N (x) + fi-l,N(x) + fi{x)T R + fi tN - i(x) - 4 fi tN {x) 

(30) 


5.1.4 Case III: | i — j\ — 1 

Finally we calculate ^f(f iti+ i(x)) for the off-diagonal elements f i)i+ i(x) = 

XiX i+ 1 - In this case 

■5?(fi,i+l (x)) = Jf L (f h2 (x)) + 2Zb(f iii+ i(x)) +Jf R (f N -i, N (x)) 

Fix i — 2, N — 2. It is easy to verify that 

2%(/,,. + i(x))=/,i + i(r i ’ i + 1 (x))+/,, i+I (r- 1 '‘(x))+/ i , i+1 (T i+1 " 2 (x))-3/,, i+I (x) 

(31) 

At the left boundary for i — 1 the generator acting on / 12 (x) is given by 
^L(fMx)) = MT^ix)) + fi 2 (T 0,1 (x)) + fi 2 (T 2,3 (x)) - 3 f 12 (x) 
and for i = N — 1 

■^R(fN-l,N(x)) = 

fN-l,N(T N 1,N (x)) + fN-l,N(T N 2,N 1 (x)) + f n~ 1 ,n(T N,N+ 1 (x)) — 3/at_i i at(x) 
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which yields 


Sf L (f lfl (x)) = {1 - X) 
-&b{fi,i+l(x)) = 
(1-A) 

■^rUn-^N^x)) = 

(1-A) 


C(fu(x) + f 22 {x)) + + T L f 2 {x)) — (1 + B)f 12 (x) 


C(fii(x ) + fi + i t i + i(x)) + + fi,i+ 2 (x)) — (1 + i(x) 


C , (/aT-1,JV-i(^) + Jnn{x)) + ~(fN-2,N(x) + Tr/jV-i(x)) — (1 + -B)/jv_i,Jv(^) 

(32) 


and 


C \= C(a, A) = -^ + a(l - A) 
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